home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dtrcon.z / dtrcon
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDTTTTRRRRCCCCOOOONNNN((((3333FFFF))))                                                          DDDDTTTTRRRRCCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DTRCON - estimate the reciprocal of the condition number of a triangular
  10.      matrix A, in either the 1-norm or the infinity-norm
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DTRCON( NORM, UPLO, DIAG, N, A, LDA, RCOND, WORK, IWORK, INFO
  14.                         )
  15.  
  16.          CHARACTER      DIAG, NORM, UPLO
  17.  
  18.          INTEGER        INFO, LDA, N
  19.  
  20.          DOUBLE         PRECISION RCOND
  21.  
  22.          INTEGER        IWORK( * )
  23.  
  24.          DOUBLE         PRECISION A( LDA, * ), WORK( * )
  25.  
  26. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  27.      DTRCON estimates the reciprocal of the condition number of a triangular
  28.      matrix A, in either the 1-norm or the infinity-norm.
  29.  
  30.      The norm of A is computed and an estimate is obtained for norm(inv(A)),
  31.      then the reciprocal of the condition number is computed as
  32.         RCOND = 1 / ( norm(A) * norm(inv(A)) ).
  33.  
  34.  
  35. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  36.      NORM    (input) CHARACTER*1
  37.              Specifies whether the 1-norm condition number or the infinity-
  38.              norm condition number is required:
  39.              = '1' or 'O':  1-norm;
  40.              = 'I':         Infinity-norm.
  41.  
  42.      UPLO    (input) CHARACTER*1
  43.              = 'U':  A is upper triangular;
  44.              = 'L':  A is lower triangular.
  45.  
  46.      DIAG    (input) CHARACTER*1
  47.              = 'N':  A is non-unit triangular;
  48.              = 'U':  A is unit triangular.
  49.  
  50.      N       (input) INTEGER
  51.              The order of the matrix A.  N >= 0.
  52.  
  53.      A       (input) DOUBLE PRECISION array, dimension (LDA,N)
  54.              The triangular matrix A.  If UPLO = 'U', the leading N-by-N upper
  55.              triangular part of the array A contains the upper triangular
  56.              matrix, and the strictly lower triangular part of A is not
  57.              referenced.  If UPLO = 'L', the leading N-by-N lower triangular
  58.              part of the array A contains the lower triangular matrix, and the
  59.              strictly upper triangular part of A is not referenced.  If DIAG =
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDTTTTRRRRCCCCOOOONNNN((((3333FFFF))))                                                          DDDDTTTTRRRRCCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              'U', the diagonal elements of A are also not referenced and are
  75.              assumed to be 1.
  76.  
  77.      LDA     (input) INTEGER
  78.              The leading dimension of the array A.  LDA >= max(1,N).
  79.  
  80.      RCOND   (output) DOUBLE PRECISION
  81.              The reciprocal of the condition number of the matrix A, computed
  82.              as RCOND = 1/(norm(A) * norm(inv(A))).
  83.  
  84.      WORK    (workspace) DOUBLE PRECISION array, dimension (3*N)
  85.  
  86.      IWORK   (workspace) INTEGER array, dimension (N)
  87.  
  88.      INFO    (output) INTEGER
  89.              = 0:  successful exit
  90.              < 0:  if INFO = -i, the i-th argument had an illegal value
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.